Scaling Approach to Calculate Critical Exponents in Anomalous Surface Roughening 
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We study surface growth models exhibiting anomalous scaling of the local surface fluctuations. 
An analytical approach to determine the local scaling exponents of continuum growth models is 
proposed. The method allows to predict when a particular growth model will have anomalous 
properties (a 7^ oti oc ) and to calculate the local exponents. Several continuum growth equations are 
examined as examples. 
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Kinetic roughening of surfaces in nonequilibrium con- 
ditions has been a subject of great interest in the past few 
years. This is mainly due to the many important applica- 
tions of the theory of surface growth including molecular- 
beam epitaxy (MBE), fluid flow in porous media and 
fracture cracks among others [Q . Growth processes have 
very often been shown to exhibit scaling properties that 
allow one to divide the models into universality classes 
characterized by the value of the critical exponents [|],|| . 
The global interface width in a system of total lateral size 
L, which is the root mean square of the fluctuations of 
the surface height, scales according to the Family- Vicsek 
ansatz [El as 



W(L,t)=t a / z f(L/t 1 / z ), 
where the scaling function f(u) behaves as 

if u < 1 



m 



ir 
const. 



if u > 1 



(1) 



(2) 



The roughness exponent a and the dynamic exponent 
z characterize the universality class of the model under 
study. The ratio j3 = a/ ' z is the time exponent. 

However, very recent numerical studies have revealed 
a rich variety of interesting phenomena, which are poorly 
understood at present. In particular, it has been shown 
that many growth models H4lll exhibit anomalous 
roughening, i.e. a different scaling for the global and the 
local surface fluctuations. This leads to the existence of 
an independent local roughness exponent ai oc that char- 
acterizes the local interface fluctuations on scales I <C L 
and differs from a. 

A general analytical method to determine anomalous 
exponents from the continuum growth equations is still 
lacking. In particular, a dynamic renormalization-group 
(RG) calculation in the Lai, Das Sarma and Villain 
(LDV) nonlinear continuum model has had only limited 
success [O] because of the existence of nonperturbative 
infrared singularities that are inaccessible to the usual dy- 
namical RG analysis. Therefore, most efforts have been 
focused on determining the anomalous critical exponents 
using simulation of discrete models or direct numerical 
solutions of the Langeving-type equations of growth. 



In this Letter, we propose a new analytical approach to 
determine the scaling exponents of the local surface fluc- 
tuations in continuum growth models exhibiting anoma- 
lous kinetic roughening. Our method allows to predict 
when a particular growth model is expected to have 
anomalous properties (a 7^ a>i oc ) and, in principle, also to 
determine the local exponents. We illustrate the method 
by studying several growth equations with and without 
anomalous scaling. 

We shall be interested here in continuum growth mod- 
els in d + 1 dimensions which dynamics is expected to be 
described by the Langevin-type equation 



dh ., ._, . . . 
— = *(Vfc)+»?(x,t), 



(3) 



where h(x, t) is the height of the interface at substrate 
position x at time t. The actual form of the function 
$(V/i) defines a particular model |13|| . r](x,t) is a noise 
term uncorrelated in space and time, (t/(x, t)r)(x', if)) = 
2D5(x-x')S(t-t'). 

The scaling properties of the local surface fluctuations 
can be investigated by computing cither the height-height 
correlation function, G(l,t) = {[h(x + 1, t) — h(x, t)] 2 ), 
where the average is calculated over all x (overline) and 
noise (brackets), or the local width, w(l,t) = (([/i(x, t) — 
C 1 )*] 2 )*) 1 ; where (• • •); denotes an average over x in win- 
dows of size I. For growth processes in which an anoma- 
lous roughening takes place these functions scale as 



w(l,t)~^G(l^j = tPf A (l/t l / z ), (4) 

with an anomalous scaling function pHy] given by 



!a{u) 



u ai °" if u < 1 
const, if m> 1 



(5) 



instead of Eq. (0) . The standard self-afflne Family- Vicsek 
scaling 0] is then recovered when a = <Xi oc - 

It is important to note that anomalous scaling stems 
from the fact that the mean square local slope has a 
non-trivial dynamics. A standard (Family- Vicsek) self- 
afflne scaling of the local interface fluctuations, i.e. a = 
aioc-, implies that the square local slope G(l = a,t) = 
([/i(x + a, t) — h(x, t)] 2 ), where a is the lattice spacing, 



becomes G(a,t) ~ constant very rapidly in time. One 
can also see that this constant must go to zero as a — > 0. 
However, as can be easily seen from Eqs.([|) and (g), in 
growth models exhibiting anomalous scaling, the local 
slopes scale as 



G(a,i)-t 2 



(G) 



where the exponent re = /3 — a/ oc / ' z > 0. The existence 
of such a diverging mean local slope introduces a new 
correlation length in the growth direction, which enters 
the scaling of the local fluctuations of the height. There- 
fore, in any growth model, the existence (or absence) of 
anomalous scaling of the local height fluctuations can be 
investigated by computing the mean local slope G(a,t). 
In the continuous limit, a — > 0, G(a,t) can be written as 
G(a,t) ~ s(t)a 2 where 



»(*) = <(v/0 5 



(7) 



This corresponds to calculating the mean square local 
derivative of the interface height. In general, this quan- 
tity scales as a power law s(t) ~ t 2K . Negative values 
of the exponent re will result in a normal Family- Vicsek 
scaling of the local fluctuations with the same rough- 
ness exponent. On the contrary, for re > the correla- 
tion length s(t) diverges in time and becomes a relevant 
length scale in the problem that changes the local scal- 
ing, as discussed above. In this case, anomalous scaling 
with a local roughness exponent 



a loc = cc — zk, 



(8) 



is expected to occur. As we will see in the following, this 
simple observation allow us to find a general method to 
compute anomalous critical exponents. Note also that 
the exponent re corresponds to the anomalous time expo- 
nent/3, inRef. ggg]. 

Now, let us see how the scaling behaviour of the mean 
local slope s(t) can be obtained from the continuum 
growth equation. By applying the operator V to both 
sides of the growth equation (||) one gets to an equation 
of motion for the local derivative T(x, t) = Vh of the 
form 



dT 5$ „„ „ 

= — VT + Vr?. 

dt ST ' 



(9) 



In general, this Langevin equation may contain nonlinear 
terms that break the translational symmetry h —* h + c 
in the growth direction. This implies that the resulting 
interface T(x, t) may not be rough (a = 0). 

The global width of the interface T(x, t) is given by 



W T (t) = (T(^t) 2 ) 1 ' 2 = {{Vh) 2 ) 1 ' 2 



(10) 



where (V/i) = has been used. This leads to the gen- 
eral result that s(t) — Wy{t) and the exponent re of the 



average local interface slope s(t) ~ t 2fi is given by the 
time exponent of the global width of the local derivative 

T(x,t), 



W T (t) ~t K 



(11) 



It then becomes clear that one could obtain the anoma- 
lous exponents of the interface h(x, t) by finding the time 
scaling behaviour of the fluctuations of T = V/i. 

In the following we investigate the existence of anoma- 
lous scaling in several continuum growth models. In the 
examples that we analyse here, a Flory-type approach in- 
troduced by Hentschel and Family [[L5[ suffices to obtain 
the exponent re from the corresponding Eq.(H) for every 
model in good agreement with existing simulation results. 
The Flory approach can be seen as a stability analysis of 
the equation of motion (ph for the corresponding local 
derivative T(x, t) of the interface. 

Kardar-Parisi- Zhang equation.- The first system we 
examine is the noise-driven interface growth model given 
by the Kardar-Parisi-Zhang (KPZ) equation |14fl 



r)h 

j^ = vV 2 h + \(Vh) 2 + v (x,t). 



(12) 



This equation was originally introduced to describe bal- 
listic deposition growth far from equilibrium. In 1 + 1 
dimensions the critical exponents a — 1/2 and z = 3/2 
can be calculated exactly 14 1. On the basis of much nu- 
merical work, it is well established that the KPZ equation 
does not exhibit anomalous roughening. Let us see how 
this result can be derived in dimension 1 + 1 by use of 
our approach. 

In this case, the growth equation for the local deriva- 
tive, Eq.(H), reads 



dT 
~dt 



d 2 T , d - dr, 

f-^r + A — T 2 + —S 



dx 2 



dx 



dx' 



(13) 



in 1 + 1 dimensions. From this growth equation the scal- 
ing behaviour of the fluctuations of the interface T(x,t) 
can be obtained. We find that the width scales as 
Wx{t) ~ i^ 1 / 5 , where the exponent can be obtained by 
a Flory-type approach |15|| as follows. 

In the spirit of Ref. |15| , we assume that at long times 
t 3> ti, and averaged over length scales /, the typical mag- 
nitude of the fluctuations in T(x,t) scale as T;, and that 
these fluctuations last for times of the order ti . The idea 
now is to estimate the magnitude of the individual terms. 
Basically for any equation such as Eq.(H3) to show time 
scaling each separate term, when coarse-grained over 
length scales I, must be of the same order of magnitude 
or negligible. Only under these circumstances can scal- 
ing behaviour arise. The various terms in Eq.(n3|) may be 
estimated as (\dT/dt\)i - T t /t h (\d 2 T/dx 2 \) l - T z /Z 2 , 
(\dT 2 /dx\)i ~ Tf/l. As for the noise, one can estimate 
its fluctuations on lenght scales I and time scales ti as 
{\d-q/dx\)i ~ (Z 3 ^) -1 / 2 for smooth surfaces, while for 



rough surfaces (\dr)/dx\)i ~ (T;/ 2 i/) -1 / 2 (see Ref. Q 
for details). Note that the correct scaling of the noise 
term depends on the particular equation under study. 

At sufficiently large length scales the nonlinear term in 
Eq.(O) will dominate the diffusion term, and so, equat- 
ing the (\dT/dt\)i term with the nonlinear term we ob- 
tain that a typical fluctuation scales as Y; ~ l/ti. To 
proceed further we now equate our estimate for the noise 
fluctuation (T;Pt;) -1 / 2 to the inertial term, which gives 
Y; ~ (ti/l 2 ) 1 ' 3 . So, we can estimate that a fluctuation 
of Y scales in time as Y; ~ t7 ' . 

So, in the case of the KPZ equation in f + 1 dimensions 
a negative exponent k = —1/5 is found indicating a stan- 
dard scaling as expected. A similar Flory-type computa- 
tion also shows that, in fact, the KPZ model exhibits a 
self-affine scaling in d+1 dimensions and k = — l/(4 + d). 

A particular case of the KPZ equation is the linear in- 
terface growth model (A = 0), first studied by Edwards 
and Wilkinson |l6| . For this model Eq. (|lj) can be solved 
exactly and we obtain the exponent k = —1/4, as corre- 
sponds to a standard Family- Vicsek scaling. 

Surface growth with conservation law.- The KPZ equa- 
tion does not conserve the total volume of the interface. 
The conserved version of KPZ was studied by Sun, Guo 
and Grant |17[] and is given by 



In this case the growth equation for the local derivative 



dh 

at 



-KV 4 h + AV 2 {Vhf + r/ c (x, t), 



(14) 



where 



(r /c (x, t)r] c (x', t')) = -2DV 2 d(x - x')<5(t - *')■ (15) 

Here the exponents are exact in any dimension fll7| ], in 
particular a — 1/3 and z = 11/3 for d — 1. We have 
investigated the possibility of anomalous scaling in this 
model in 1 + 1 dimensions. From the corresponding equa- 
tion for the local derivative Eq.(||) and the noise fluc- 
tuations ~ (Z 3 £i) -1 / 2 , we obtain the scaling behaviour 
W r (t) ~ t~ 2/u for the fluctuations of T(x, t). This result 
shows that the SGG equation has also a normal scaling 
of the local fluctuations of the height. 

Linear MBE model.- As a simple example of anoma- 
lous roughening, we now consider the linear model for 
MBE growth pUlJJ given by 



dh 



-KV 4 h + n{x,t), 



(16) 



Despite its simplicity, this equation has played an impor- 
tant role in the theory of MBE. The critical exponents 
are easily calculated in any dimension, and in particu- 
lar one has a = 3/2 and z = 4 in dimension d = 1. 
The model exhibits super-rough interfaces (a > 1). This 
leads to anomalous (super-rough) scaling flqfl. In fact, 
numerical simulations pq , p9| ,pp|] of Eq. ( jlq) showed that 
the local scaling is given by Eqs.fl) and (ra) with a local 
roughness exponent a>i oc ~ 1. 



<9Y d A T d , , 

-dt = ~ K W + d^ r]{x ' th 



(17) 



is linear and can be easily solved. The Flory approach 
now gives the exact exponent n. We can estimate the 



curvature diffusion term as {\d T / dx 



Ti/F, being 



the estimate for the noise term {\dn/dx\)i ~ (Z 3 ^)^ 1 / 2 
in this case. Equating fluctuations of each of the two 
terms on the right-hand side of Eq. ( |17| ) with the inertial 
term, we obtain the scaling behaviour Wy(t) ~ t 1 ^ 8 for 
the width of the local interface derivative. The positive 
value of the exponent k means that the local slope s(t) 
becomes a relevant correlation length in the problem and 
anomalous roughening is to be expected. The scaling re- 
lation (||) gives us an exact determination of the local 
roughness exponent ai oc = 1 for this model in dimension 
d = 1 . This is good in agreement with existing numerical 
results |l|Ul|,|]. 

Random diffusion model.- So far we have been consid- 
ering growth models in which the only source of random- 
ness is in the influx of particles on the surface. Let us 
now consider the growth model 



dh d r _, . dh, , ,. 

m=d- X [D{x) d-x ]+ ^ t) > 



(18) 



where D(x) > is a quenched columnar disorder with 
no correlations. This random diffusion coefficient is dis- 
tributed according to P(D) ~ D~^ for D < 1 and 
P(D) = for D > 1. This model describes a fluid inter- 
face advancing through a stratified porous medium with 
long range correlations in the growth direction. 

The random diffusion model was originally introduced 
|6j to demonstrate that anomalous roughening is not due 
to super-roughening effects, but it can also take place in 
models in which a < 1. The critical exponents depend 
on the disorder strength fi and can be calculated exactly 
§ a = l/[2(l-/i)] and z = (2-/x)/(l-/i) for < n < 1. 
This model exhibits anomalous roughening with a local 
roughness exponent a/ oc = 1/2 that is independent of the 
disorder |^|. 

The existence of anomalous roughening in this model 
can be rederived by use of our approach as follows. The 
growth equation for the local derivative 



<9Y d 2 d 

m=^ D ^ + d- x ^ 



(19) 



can be solved exactly to obtain the time scaling behaviour 
of W T {t) ~ t K . We find that k = /i/[2(2 - n)] for 
< /i < 1, which gives us a local roughness exponent 
etioc = a — zk = 1/2 in agreement with earlier results [pi. 
Lai- Das Sarma- Villain equation.- The last example we 
study is the LDV @|IJ equation for MBE growth 



dh 
~dt 



-KV 4 h + AV 2 {Vhf + 77(x, f). 



(20) 



This equation is expected to describe the behaviour of 
the long-wavelength fluctuations of the interface height 
in several discrete models of MBE growth (HQ^^Q. 
Note that this equation differs from the SGG equation 
discussed above in the non-conserved character of the 
noise in this case. 

According to a dynamical RG analysis of Eq. (E(l) the 
global exponents a = (4 — d)/3 and z = (8 + d)/3 are ex- 
pected to be exact to all loops. However, numerical simu- 
lations of the LDV equation and several discrete growth 
models in the same universality class have shown that 
the model exhibits anomalous scaling with a local rough- 
ness exponent ai oc ~ 0.7 P,[7Ull|l in 1 + 1 dimensions. 
A perturbative dynamic RG approach to Eq. (|20| ) showed 
|l2| the existence of a strong-coupling infrared singular- 
ity that prevented from obtaining the local anomalous 
behaviour by perturbative methods. 

Again, the existence of anomalous scaling in this model 
in 1 + f dimensions can be investigated by use of the 
growth equation for the local derivative of the interface. 
In this case Eq.(g) reads 

<9T 
~dt 



d 4 T d 3 

dx A dx 3 



T 



9 , , 



(21) 



A Flory approach can also be applied to this case 
to determine the scaling of Wr ~ t K . As we did 
for the KPZ equation, we can estimate the terms in 
Eq.© as <|0T/0t|>, - Ti/ti, (\d*T/dx% ~ T,/Z 4 , 
(\d 3 T 2 /dx 3 \)i ~ Tf/l 3 , and the noise being (\drj/dx\)i ~ 
(T/Z 2 £i) -1 / 2 as in the KPZ case. Assuming that the non- 
linear term dominates the curvature diffusion term at 
large scales, we obtain Wy ~ t 1 / 11 . 

According to this the exponent k — 1/11 is positive 
and, as a consequence, the LDV equation in 1 + 1 dimen- 
sions is expected to display anomalous roughening. The 
local roughness exponent ai oc = 8/11 ~ 0.73 is given by 
the scaling relation (0) with a = 1 and z — 3 for d = 1. 
Our estimate for a/ oc is in excellent agreement with the 
numerical result a/ oc ~ 0.73 ± 0.04 ||. 

Conclusion.- In summary, we have introduced a new 
method to calculate the local roughness exponent in 
growth models exhibiting anomalous kinetic roughening. 
We have shown that a divergent dynamics of the mean 
square local derivative, s(t) ~ t 2K with k > 0, is the 
cause for an anomalous roughening of the local surface 
fluctuations. The exponent k can be obtained by study- 
ing the time scaling behaviour of the fluctuations of the 
local derivative. The local roughness exponent can then 
be obtained by use of a scaling relation, Eq.(||). We have 
examined the existence of anomalous roughening in sev- 
eral models. For the linear models our results are exact. 
Nonlinear growth equations were studied by a scaling ap- 
proach, which gave results in good agreement with exist- 
ing simulations. 

Finally, our findings indicate that a self-affine (a — 
ai oc ) Family- Vicsek scaling of the interface fluctuations 



is the result of a fine balance among the relevant terms in 
the growth equation. This balance can be altered by in- 
cluding conservation laws or higher order nonlinear terms 
resulting in a divergent behaviour of the local derivatives 
s(t) and anomalous scaling. 
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